Vibrated polar disks: spontaneous motion, binary collisions, and collective dynamics 
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We study the spontaneous motion, binary collisions, and collective dynamics of "polar disks", i.e. 
purpose-built particles which, when vibrated between two horizontal plates, move coherently along 
a direction strongly correlated to their intrinsic polarity. The motion of our particles, although 
nominally three-dimensional and complicated, is well accounted for by a two-dimensional persistent 
random walk. Their binary collisions are spatiotemporally extended events during which multiple 
actual collisions happen, yielding a weak average effective alignment. We show that this well- 
controlled, "dry active matter" system can display collective motion with orient at ionally-ordered 
regions of the order of the system size. We provide evidence of strong number density in the most 
ordered regimes observed. These results are discussed in the light of the limitations of our system, 
notably those due to the inevitable presence of walls. 

PACS numbers: 64.60. Cn, 05.70.Ln, 05.65.+b, 45.70. Qj 



I. INTRODUCTION 

Active matter is nowadays the expression designating 
out-of-equilibrium systems where energy is spent locally 
to produce directed motion pQ. Examples abound, at all 
scales and in many fields, with perhaps the most spectac- 
ular ones in living systems, from the motion of subcellular 
components to that of large animal groups [2HTU]. Ac- 
tive matter is of course the word used by physicists, who 
have been flocking to this rather new, burgeoning field, 
producing novel theoretical ideas and results, and creat- 
ing well-controlled objects that can move "on their own" 
either swimming in a fluid or crawling on a substrate [TT1- 
[T7] . Although the fluid is always there, it can sometimes 
be safely neglected, and momentum conservation thus 
disobeyed, in particular when objects are moving on a 
substrate, which acts as a momentum sinks. One then 
speaks of "dry" active matter [TJ [18] . 

Microscopic models of dry active matter usually con- 
sist of self-propelled particles whose local interactions 
give rise to some kind of alignment. This includes the 
now famous Vicsek et al. model in which constant-speed 
point particles align "ferromagnetically" with neighbor- 
ing ones, in the presence of noise, in what amounts to an 
out-of-equilibrium version of the XY model [19]. Impor- 
tant results have been obtained recently both on micro- 
scopic models like the Vicsek model and on continuous 
descriptions of dry active matter [20-26 . These include 
the seminal calculation by Toner and Tu [27] [28] who 
confirmed, via a renormalization-group approach, the nu- 
merical findings of Vicsek et al. that collective motion 



can be sustained at all scales in spite of the presence of 
noise: true long-range orient at ional order exists even in 
two dimensions when spins are moving in space. 

Although much remains to be done, progress has been 
recorded in matching the many numerical results ob- 
tained with microscopic models to their continuous de- 
scriptions [24H261 [29] , and important theoretical predic- 
tions have been confirmed in numerical studies [2TH23] . 
The situation is much less satisfactory at the experimen- 
tal level. This should not be a surprise: living systems, 
be they large animal groups or collections of cells or even 
biofilaments and motor proteins, are difficult if not im- 
possible to control, and they typically involve many ex- 
tra, unwanted, and often unknown factors and mecha- 
nisms. 

In this setting, building inert objects capable of mov- 
ing by themselves or by capturing energy from their en- 
vironment is a crucial step towards well-controlled exper- 
iments. In the context of dry active matter, two options 
are available: moving "robots" interacting by direct con- 
tact or at some distance via sensors [30], or asymmetric 
objects vibrated on some homogeneous substrate. Here 
we follow this last option, building on knowledge gath- 
ered over the last few decades on shaken granular media. 

Yamada, Hondou, and Sano were pioneers in demon- 
strating that an axisymmetric polar object vibrated be- 
tween two plates can move quasiballistically [31]. Ku- 
drolli's group studied the behavior of polar rods [32 and, 
more recently, of short snakelike chains [33] . Shaken elon- 
gated apolar particles (a realization of so-called active 
nematics) have been considered in [34] . 
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FIG. 1. (color online) Vibrational system, (a) The vibration is provided by an electromagnetic servo-controlled shaker ballasted 
with lead and isolated from ground by rubber mats (4). It is then transmitted to the vibrating plate (2) via a brass rod coupled 
to a square air bearing (3). The vibrating plate (2) is composed of a nylon plate, an expanded polystyrene truncated cone, 
on top of which is glued the glass plate which supports the vibrated grains. Finally, the grains are confined vertically by a 
top glass plate fixed to an external rigid frame (1). (b) This top glass plate (1) is in a rigid frame (2), the vertical position of 
which is finely tuned by micro-control screws (3) to ensure an homogeneous gap of 2.4 mm. (c) The vibration is imposed and 
retro-controlled at the center of the glass plate supporting the grains. It is measured at the periphery of the plate. Vertical 
and azimuthal transfer functions between the acceleration imposed by the controller in the center of the plate and the one 
measured at six points on the plate perimeter. The horizontal dashed dark lines indicate a control of the acceleration within 
5%. The vertical dot-dashed line indicates the frequency at which all experiments reported in this paper have been conducted 



The particles above are all elongated objects, with 
their asymmetry inscribed in their shape. This usually 
ensures that their (inelastic) collisions result in align- 
ment, in a way after all rather similar to what happens in 
numerical models [35H38] . In this paper, we deliberately 
depart from this: we report further and at length on our 
experiments performed on vibrated polar disks, follow- 
ing a first account of their collective properties [39 . The 
asymmetry of our particles is thus not in their shape, but 
in their contact properties with the vibrating plate (both 
friction and restitution coefficient), which endows them, 
as we shall see, with rather unusual collision dynamics. 

The rest of this paper is organized as follows: in Sec- 
tion [n] we detail our vibration apparatus and the design 



tion of the individual motion of our particles and how 
it varies with our main control parameter, the vibration 
amplitude. In Section |IV[ we show that collisions be- 
tween two of our polar disks possess rather complicated 
properties, but they do result in some effective alignment. 
In Section |V| we turn to the collective properties of our 
particles and show that our system can produce orienta- 



tional order over large scales. Section [VI] contains a dis- 
cussion of the limitations of our results and offers some 
perspectives on future work. 



of our polar disks. Section III is devoted to the descrip- 



II. EXPERIMENTAL APPARATUS 

In principle, producing some kind of self-propulsion us- 
ing vibrating man-made objects allows for incomparably 
more control than with, say, biological organisms or even 
motility assays in which biofilaments are displaced by 
motor proteins grafted to a substrate. Nevertheless, in 
order to be free from various artifacts, these experiments 
require extremely homogeneous and strictly horizontal 
vibration, which is notoriously difficult to achieve. In- 
deed, previous attempts have been suspected of being 
plagued by large-scale spurious modes [3H EQl SB- In 
this section, we first introduce the vibrational system we 
used to generate a well-controlled vertical vibration. We 
then describe and motivate the design of our polar disks. 



A. Vibration system 
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FIG. 2. (color online) Self-propelled polar disks, (a) Side 
and bottom views of a polar disk with the built-in polarity 
n. The white part of the particle is made of copper-berylium, 
while the grey part is made of nitrile. (b) Side and top views 
of the polar disks with their respective polarities. The black 
scale bar is 4 mm. 



The vibration is provided by an electromagnetic servo- 
controlled shaker (V455/6-PA1000L,LDS) [FigjTJi.], 
which rests on a thick wooden plate ballasted with 300 kg 
lead bricks and isolated from ground with rubber mats 
(MUSTshock 100xl00xEP5). The vibration is transmit- 
ted to the vibrating plate via a 400 mm- long, 8 mm-thick 
brad ross, coupled with a stiff square air-bearing slider of 
size 127 mm (C40-03100-100254, IBSPE). The brass rod 
is flexible enough to compensate for the alignment mis- 
match, but stiff enough to ensure the mechanical coupling 
between the slider and the shaker. The slider ensures vir- 
tually friction-free and submicron amplitude horizontal 
motion. The plate is composed of a 50 mm thick nylon 
plate fixed to the slide, a 50 mm expanded polystyrene 
truncated cone, on top of which is glued the 425 mm di- 
ameter glass plate supporting the vibrated grains. Fi- 
nally, the vertical motion of the grains is confined by 
a top glass plate attached to an external non vibrating 
rigid frame. The two glass plates are separated by a gap 
h = 2.4 mm. 

An uniaxial accelerometer, used to control the vibra- 
tion, is glued at the bottom of the bottom glass plate, 
inside a small cavity encarved in the polystyrene. In or- 
der to control the vibration quality, we measure with a 
triaxial accelerometer (356B18, PCB electronics) the ver- 
tical and azimuthal transfer functions between the con- 
troller and the acceleration recorded on the perimeter of 
the bottom glass plate, when vibrated at an acceleration 
Aq = 1 g, where g stands for the acceleration of gravity 
[FigJTJ).]. A working frequency / = 115 Hz was chosen 
so that it avoids all resonances in the vibrating system 
while inducing motion of our polar disks (see below). At 
this frequency, the azimuthal to vertical acceleration ra- 
tio is smaller than 1% and the relative vibration hetero- 
geneities are less than 5%. 



B. Self-propelled particles 



To the best of our knowledge, all previous experiments 
looking at the coherent displacement of objects collid- 
ing with a vertically-vibrated horizontally-homogeneous 
substrate have used elongated particles [3TJ [32] . In other 
words, the axis necessary to break rotational symmetry 
was simply encoded in the shape. In such a case, col- 
lisions directly induce alignment. Here, by contrast, we 
use essentially circular particles: they consist of a top 
metallic disk with two asymmetric "legs" aligned to give 
the object a polar axis (Fig. [2^i). These particles inter- 
act strictly via their top circular hard parts, much like 
more conventional hard disks. Their polarity lies in the 
contacts with the vibrating plate. 

In detail, we micro-machined copper-beryllium disks of 
diameter 4 mm and height 2 mm with an off-center leg 
to which a nitrile skate of hardness 90 Sh was glued at 
a diametrically opposite position (Fig. . The skate is 
overmoulded to the copper-beryllium thanks to two poly- 
meric layers, one adhesive to the metal (THIXON(TM) 
520 Adhesive, Rhom and Haas), the other to the ni- 
trile (THIXON P-11/0285SPL/25KG, Rhom and Haas). 
These two "legs", which have different mechanical re- 
sponse under vibration, endow the particles with a polar 
axis which can be determined from above thanks to a 
black spot located on their top (Fig. ^p). 

In the following we also use plain disks (same metal, 
diameter and height) as a control to ensure that all 
the effects reported here are indeed due to the "self- 
propulsion" induced by the built-in polarity. 
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C. Data acquisition 

The position and polarity of the disks were recorded 
by a CCD camera with a spatial resolution of 1728x1728 
pixels. The acquisition rate was set to 20 Hz or 25 Hz, 
i.e. about five times slower than the shaking frequency of 
115 Hz. The ratio of the vibration frequency and the data 
acquisition rate sets the value of tq, the smallest time in- 
crement accessible. The trajectory of each particle was 
then reconstructed using standard tracking software. In 
the following, the time unit is set to be the vibration pe- 
riod and the unit length is the particle diameter. Within 
these units, the resolution on the position r of the parti- 
cles is better than 0.05, and that on the orientation ft is 
of the order of 0.05 rad. 



III. SINGLE PARTICLE MOTION 

Our polar disks, once vibrated at / = 115 Hz in our ap- 
paratus, undergo fairly coherent, persistent motion of the 
particles, over some range of easily accessible amplitudes 
T = A uj 2 /g with u = 2nf. 



A. Typical behavior 

The single particle motion, which arises from the inter- 
play of collision/friction of the two legs with the bottom 
plate, and of the metallic disk with the top plate, is too 
complicated to be described quantitatively. In the follow- 
ing, we only discuss its projection on the horizontal plane, 
"forgetting" about the three-dimensional nature of the 
events at its origin. Typical two-dimensional trajectories 
of isolated particles appear smooth and rather straight 
(Fig ,3i), as opposed to those of the control plain disks 
(Fig. 3d). A closer inspection, though, reveals that our 
polar particles typically move steadily forward for some 
time, but also intermittently stop or even move backward 
briefly (Fig {3^) and Supplementary movie [42]). 

The directed, persistent motion of particles is recorded 
in the (discrete-)time series of their position r(t) and po- 
larity 0(t). The polarity is very persistent in time: the 
distribution of A0(r) = Q(t + r) — 6(t), its increment over 
time r, is essentially exponential and sharply peaked at 
zero for t = tq. Moreover, the distribution of the reduced 
increment A0(r) = A^(r)/(A^(r) 2 ) 1 / 2 remains roughly 
independent of r up to fairly large values of r, after which 
it gradually recovers the expected asymptotic Gaussian 
character (Fig.|3]i). The crossover value of r (of the order 
of 100) can be taken as a rough estimate of the persis- 
tence time of polarity. 

By and large, the displacement of the particles is over- 
whelmingly along their polarity axis: the distribution 
function of the angle a(r) between Ar(£, r) = f(t + r) — 
r(i), the displacement vector over time r, and 0(t) is 
peaked around zero. Strikingly, this is all the more so 
as r is large: The backward motion events apparent in 
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FIG. 3. (color online) Statistical properties of the individ- 
ual motion for r = 2.8. (a) Sample trajectories of the self- 
propelled particles. The red arrows indicate the instantaneous 
orientation of the polarity. The black arrows indicate the ori- 
entation of the displacement between two successive frames, 
(b) Sample trajectories of the isotropic particles. The black 
arrows indicate the orientation of the displacement between 
two successive frames, (c) Five time series of the displacement 
component along the polarity. The signals have been shifted 
for clarity. The dotted lines indicate the zero. One clearly 
observes negative events, corresponding to backward motion, 
(d) Distribution of the reduced increments of the polarity ori- 
entation A0(t) over a time r for increasing r. (e) Distribution 
of the angle a(r) between the displacement vector Ar (£, r) 
and the polarity n(t) for increasing r. (f) Distribution of 
| Ar(£, t)|/t, the displacement over a time r, normalized by r 
for increasing r. The color code in frames (d,e,f) is indicated 
by the legend in frame (f). 



the secondary peaks at ±tt present at small r average 
out as r is taken larger, an indication that during most 
of these "backward events" the polarity remains largely 
unchanged (Fig.pe). Again, this is only true up to some 
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crossover value of r ~ 100 beyond which the distribution 
of a(r) must gradually become flat. For r larger than the 
crossover, one indeed observes a widening of the distri- 
bution. This displacement mostly along the polarity axis 
is performed at a fairly well-defined speed: for not too 
large r, the distribution of |Ar(t, r)\/r is peaked around 
a most probable value. Increasing r from To to about the 
crossover time mentioned above, the distribution keeps 
the same most probable value and gets narrower and nar- 
rower (Fig. |3]F). The well-defined most probable value is 
thus nothing but the average speed (v). This indicates 
again that over these timescales our particles essentially 
go straight. The distribution is then essentially Gaus- 
sian. For timescales larger than the crossover, there is 
of course a shift of the "speed" towards lower values as 
expected when the particles enter the long-time, uncor- 
rected, diffusive regime. 



B. Influence of the vibration amplitude F 

We have seen above that our polar disks can be faith- 
fully described, over scales which average out the stop- 
ping and backward events, as moving at a well-defined 
finite speed (v) while being subjected to weak rotational 
diffusion. A direct and accurate measure of the (rather 
long) persistence length/time of the trajectories of our 
particles via, say, the time decay of the autocorrelation of 
their polarity, is rendered difficult by the relatively small 
size of our system. The data presented in Fig. [3] were ob- 
tained in a dish of diameter 40 vibrated at an amplitude 
T = 2.8. In such conditions, the rather straight trajecto- 
ries will hit the wall long before their have turned enough 
to yield a significant decay of the polarity autocorrela- 
tion. To overcome this difficulty, we used the distribution 
of (normalized) polarity increments as shown in Fig. [3]i. 
Being independent of r at small r, and the mean square 
angular increment being linear in r (Fig. |4^i), it allows 
to define the rotational diffusion constant Dq as its (half- 
variance, and thus the persistence length as £ = (v) /Dq. 

The influence of Y on V m , Dq, and £ is described in 
Fig. [4] First of all, we observe a rather sudden drop in 
self-propulsion when T is decreased (Fig. [4^): for small 
amplitude, the drive is too weak and our polar disks 
do not move much. Over the range T G [2.7,3.8], on 
the other hand, the average speed varies little. In con- 
trast, the rotational diffusion constant Dq, which is easily 
captured from the root mean square angular increment 
computed on a lag r (Fig. |4Jd), steadily increases with 
T, with an almost linear behavior in the T-range over 
which steady propulsion occurs (Fig. [4^). This indicates 
that the main effect of T is on the strength of the an- 
gular noise. As a consequence, the persistence length £ 
becomes large as T is decreased within the steady propul- 
sion range [2.7, 3.8], but falls back for too small T values 
(Fig#). 
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FIG. 4. (color online) Self propulsion (a) Characteristic speed 
(v) defined by the most probable value of \\Ar(r)/r vs. T for 
r < 100. (b) Mean square angular increment computed on a 
lag time r versus r, for vibration amplitude T G [2.7 — 3.7]. 
(c) Diffusion coefficient Dq of the orientation of the polarity 
vs. r. (d) Persistence length £ = (v)/Dq vs. Y. Symbols (+) 
and (o) indicate two different experimental runs with differ- 
ent experimental conditions, such as ambient humidity, fine 
positioning of the top glass plate, etc.... The data analyzed in 
this paper for the small system correspond to the run labelled 

(+)- 



IV. BINARY COLLISIONS 



Having characterized the motion of our particles, we 
now turn to a statistical description of their binary col- 
lisions using trajectories recorded when only few parti- 
cles are evolving in the system (typically 50 particles in 
a domain of diameter ~ 40). As already stated, most 
models for the collective motion of objects moving on a 
substrate involve more or less explicit alignment rules. 
For elongated, rod-like particles, inelastic collisions im- 
mediately lead to alignment. Here, given that our par- 
ticles interact by collisions of their circular top part, no 
such direct alignment occurs. Rather, encounters typi- 
cally consist of multiple collisions each followed by some 
rebounds (Fig. |5^i and Supplementary movie [43]). It is 
thus by no means clear a priori that even some effective, 
average alignment takes place during collisions of two of 
our particles. 
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FIG. 5. (Color online) Collisions, (a) Two samples of encounters exhibiting distinct effective alignment properties. In the 
case of the encounter on top, the two particles encounter face to face, collide successively, while their polarities align, until the 
encounter ends at a relatively large distance from the first collision position. In the case of the bottom one, the two particles 
again encounter face to face. However the particles polarities essentially do not change during the successive collisions and he 
encounter ends when the impact factor of the collision reduces to zero. At the end of the encounter the particles escape back 
to back, (b) Sketch of an encounter with the definition of its geometrical properties, (c) Mean duration of the encounters as 
a function of the distance d c used as a criterion to define them, (color code as indicated in the legend); the vertical dash line 
indicate the criterion we have kept to define encounters, (d) Distribution of encounter durations r c , (color code as indicated in 
the legend), (e) Distribution of encounter lengths Ar c , (color code as indicated in the legend), (f) Average encounter duration 
< r c >, (□) and length < Ar c >, (o) as a function of I\ 



A. Defining collisions, their duration, and their 
spatial extent 

Given that encounters of two of our particles typi- 
cally involve many actual collisions, sometimes occuring 
more frequently than our sampling rate, we studied the 
statistics of encounters defined by the total space-time 
event during which two particles, initially far apart, have 
had their centers stay closer from each other than some 
threshold distance d c . To insure that these events are 
"real" collisions, this distance criterion was supplemented 
by the condition that the polarities of the two particles 
at the beginning of the encounter indicate that their dis- 
tance is likely to decrease. 

Obviously, d c must be chosen larger than 1, but not 
too large otherwise encounters may not cease. To use 
a unique, objective value of d c , we calculated (r c ), the 
mean duration of encounters, as a function of d c , over a 
set of a few thousands events. The resulting curve shows 
a plateau-like behavior in the range d c G [1.5 — 1.9], the 



middle of which can be defined as being the inflection 
point d* c± 1.7 (Fig. [5^). Furthermore, this well-defined 
value, around which the statistics of encounters do not 
vary significantly, does not vary much when T is changed 
within the useful range [2.7, 3.8]. It is thus adopted in the 
following to define encounters (hereafter renamed colli- 
sions for simplicity) quantitatively. Once collisions thus 
defined, their spatial extension Ar c can be estimated 
as the distance between the mid point of the particle 
centers at the beginning and at the end of the collision 
(Fig. ^jp) We find the collision sizes distributed roughly 
algebraically with an exponent of the order of 3, whereas 
the collision durations r c are distributed roughly expo- 
nentially. The influence of T on these distributions re- 
mains rather weak, with experimentally-observed means 
decreasing slowly with T (Fig. pHi-f). 
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FIG. 6. (color online) Alignment properties of the collisions, 
(a) Distribution of the deviation angle /3 of the collision ori- 
entation, (see Fig. [5]d) (b) Sketch of the collisions in the plane 
(2\0 in \,20 ou t) (c) Distribution of 2\Q out - in \ (d) Scatter plot 
of (2\0i n \, 20 out) characterizing the alignment properties of the 
collisions. 



or "nematic" interactions similar to that represented in 
Fig. [5|i (bottom) can be seen, but they remain rare. 

Thus, all in all, about 70% of binary collisions amount 
to no effective interaction. Among the 30% remaining 
ones, most are strongly aligning, with a few outliers. So 
far, our system seems to be rather close to the Vicsek 
model: our particles can be considered, on some coarse- 
grained timescale, to be moving at constant speed. Via 
the shaking amplitude T, we have good control on the 
effective rotational diffusion constant Dq. And binary 
collisions, when they significantly alter the orientation of 
particles, are clearly aligning. 

This said, it is still by no means clear, in the higher 
density regimes where collective motion could arise, that 
any of these observations is really important. Indeed, one 
expects there short mean free paths, and very compli- 
cated, possibly unending collision events involving more 
than two particles, a usual feature of dense collections of 
finite-size granular particles made even more acute here 
by the fact that binary collisions are spatio-temporally 
extended events. 



V. COLLECTIVE MOTION 

We now turn to the study of the collective dynamics 
of our vibrated polar disks. We first discuss interactions 
between our particles and the lateral walls of our appa- 
ratus, as well as the choice of the size and shape of the 
horizontal domain in which they are evolving. 



B. Alignment properties 

In spite of their complicated space-time structure, and 
of the intrinsically chaotic/noisy character of the parti- 
cle dynamics, the binary collisions possess a remarkable 
property: the sum of the polarities of the two particles 
changes very little between the beginning and the end 
of the collision. The distribution of the angle (3 char- 
acterizing this (see Fig. [5Jd) is sharply peaked around 
(Fig.§i). 

In the following, we thus assume that /3 = 0, and ne- 
glect the dependence of the collision statistics on the im- 
pact coefficient, so that all the information related to 
alignment is encoded in the conditional probability of 
#out, the outgoing angle of the polarities, on #i n , the in- 
coming angle (Fig. ^p). All this information is repre- 
sented in the scatter plot of out vs 6- in (Fig.[6]3,d). A first 
striking fact is that most collisions actually do not change 
polarities at all: most points of the scatter plot are lo- 
cated near the diagonal # out = — 0- in ] the distribution of 
|#out — #in| is sharply peaked near zero (Fig. [6^). In the 
(#in 5 #out) scatter plot, apart from the highly-populated 
diagonal region, most of the remaining points are located 
near the 6 out = axis (near-perfect, Vicsek- style align- 
ment) or in the triangular region between this axis and 
the 

$out — — $in diagonal (reduced outgoing angle, but 
not perfect alignment) (see Fig. |6]3,d). Ant i- alignment 



A. Domain size and shape, collisions with 
boundary 

When one of our particles hits a wall head-on, it 
bounces on it repeatedly, gradually turning after each 
hit, and eventually escapes, typically along the wall (see 
a sketch in Fig.[7£i). Thus, in a many-particle system, one 
expects an accumulation of particles near the wall, be it 
flat or circular. This is a common problem when dealing 
with self-propelled objects, as already experienced, e.g., 
by Kudrolli et al. [32]. In the "real world", one cannot 
easily rely on the periodic boundary conditions available 
in numerical studies, short of having particles move on 
some closed surface like a sphere or a torus. 

In the present case, we designed "flower- shaped" do- 
mains (Fig. ^}pc) to help particles escape from the wall. 
Crawling along a "petal", they are likely to be reinjected 
in the bulk when reaching the cusp point marking the 
limit with the next petal: each arc behaves as a launch- 
pad. Two different flower-shaped arenas have been used 
in the following, one of internal diameter 40 and 8 petals 
(the "small" system) the other one of internal diameter 
90 and 16 petals (the "large" system). The "petal" size 
was not scaled with the internal diameter, but was kept 
roughly constant: their arc radii are 9.3 for the small 
system, and 10.125 for the large one. Therefore, only the 
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FIG. 7. (color online) Boundary conditions: (a) Sketch of the 
mechanism responsible for the accumulation of self-propelled 
disks at the walls, (b) Small system: internal radius (dashed 
circle) is 20 and the arcs radii are 9.3. (c) Large system: 
internal radius (dashed circle) is 45 mm and the arcs radii are 
10.125. 

general available area increases when changing from the 
small to the large system, not the length-scale used to 
re-inject particles toward the center. 



B. Packing fraction and region of interest 

In the vast majority of models for collective motion, 
the ordered, collectively-moving phase is reached by ei- 
ther increasing the persistence length of the single parti- 
cle walk (e.g. by reducing the angular noise strength) or 
increasing the density of particles (which make interac- 
tions frequent enough to sustain the memory of the last 
aligning collision until the next event). 

In our system, we have seen that the persistence length 
is controlled rather nicely by the vibration amplitude T. 
Clearly, this is the easiest control parameter to use here, 
since changing density would be tedious. But we already 
showed that the polar disks in our system can only move 
coherently on a rather restricted range of T values (in 
particular propulsion deteriorates quickly when T < 2.7). 
This limit, in turn, imposes that we work at the "right" 
density of particles so that the transition to collective 
motion and the ordered phase be observable. If the den- 
sity is too low, our system cannot order, even for T = 2.7. 
If it is too high, on the other hand, one is bound to run 
into the problem of glassy dynamics, crystallization and 
jamming: because of their finite size, our polar particles 
flock in packed clusters. Jamming of active particles is an 
interesting problem (on which research has just started 
[441 S3) 5 but it is beyond the "simpler", Vicsek-style 
problem of collective motion of point particles considered 
here. 

Thus, that collective motion is observable in our sys- 
tem is not guaranteed. We found, by trial and error, 
that nominal packing fractions (j) (the total area covered 
by the particles in the system divided by the domain 
area) in the range <j) G [0.4, 0.6] give the best results. For 



instance, 890 particles in our "small" domain (0 = 0.47) 
vibrated at Y = 2.7 move collectively on large scales in a 
spectacular way (Fig. [8^i and Supplementary movie [46]). 

In fact, the issues discussed above are further compli- 
cated by the effect of boundaries. Even if our flower- 
shaped domains help solve the problem of particle accu- 
mulation along the walls, they do not prevent it entirely, 
and we do observe, on average, higher densities near the 
wall (Fig. 18b) . This forces us to define a "region of inter- 
est" (ROI) over which the observed time-averaged den- 
sity remains approximately constant. From graphs like 
that of Fig. [8)3, we see that the time-averaged density 
is constant to a few percent if the ROI is chosen to be 
of diameter 20, i.e. half the nominal diameter of our 
small system. For the range of packing fractions we will 
discuss below, we could safely use the same ROI size in- 
dependently of the nominal packing fraction. Finally, 
going from our small domain (diameter 40) to the large 
one (diameter 90), we find that the ROI roughly doubles 
(from diameter 20 to 40) if one keeps the requirement 
of a few percent variation of the observed time-averaged 
density. Inside the ROI, the well-defined averaged pack- 
ing fraction (/>roi is not only significantly smaller than 
the nominal packing fraction 0, but it decreases slowly 
with T (see inset of Fig|8]i). (Note that these are op- 
posite effects: decreasing T should increase order since 
this reduces Dq, but lower (/>roi values should decrease 
it.) For the 890 particles in our small system, we observe 
that 0roi e [0.36,0.42] for T e [2.7,3.7] (Fig.^). 

With these remarks in mind, we now present results 
from three sets of experiments performed by varying T 
at fixed "geometrical conditions" (Table |l|. 



TABLE I. Main characteristics of the sets of experiments 
described. 



Experiment 


Arena used 


# of particles 





typical 0roi 


Ea 


small 


890 


0.47 


0.39 


Eb 


large 


3830 


0.42 


0.26 


Ec 


small 


1090 


0.58 


0.19 



C. Onset of collective motion 

In experiments Ea and Eb, decreasing the vibration 
amplitude Y brings the system from a disordered state to 
regimes where particles, in large regions of the system of 
the order of the ROI diameter, move in approximately the 
same direction (Fig. [8^1 and Supplementary movie [46]). 

Different order parameters can be denned to quantify 
the degree of orient at ional order. The average polarity 
is easy to measure, but we preferred to use the average 
orientation of the displacement of each particle over a 
time r = 50, taken as a proxy for velocity. Although 
both choices yield essentially the same results, the latter 
one is more directly related to the transition to collective 
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FIG. 8. (color online) Emergence of collective motion: (a) Snapshot taken in a regime with long-range collective motion. The 
color code is red (reps, blue) for perfectly aligned (reps, anti-aligned) neighbors, (b) Time-averaged packing fraction (0(r, t)) t 
versus the distance r to the center of the vibrating plate for T £ [2.7 — 3.7]. The size of the region of interest (ROI) is chosen 
in order to ensure uniformity of the packing fraction inside the ROI and indicated by the vertical dashed line. The ROI in the 
small system has a radius of 10 particles diameters. We also checked that the packing fraction inside the ROI is stationary. 
(Same color code as in (c)) (c) Distribution of the order parameter ^ for T G [2.7 — 3.7], (color code as indicated in the legend), 
(d) Most probable value of the order parameter as a function of F for both the smaller experiment Ea (black o) and the 
larger one, Eb, (red +). Inset: $roi as a function of T for the same experiments, (e) Temporal evolution of the order parameter 
ty(t) for four different vibration amplitudes F = 3.6, 3.3, 3.1, 2.8. 



motion. Below, ^(t) in fact represents the modulus of 
this velocity-based order parameter, properly normalized 
to be 1 for perfectly aligned particles. 

Time series of \I/(£) show that decreasing F from 3.7 
to 2.7, the system reaches more and more ordered states, 
albeit never perfectly ordered ones as expected in a fi- 
nite system with hard walls (Fig. [8^). The most or- 
dered regimes exhibit very large fluctuations of the or- 
der parameter over very long timescales as shown also 
in the corresponding histograms of values taken by ^(t) 
(Fig#). 

As a matter of fact, our limited range of usable F values 
did not allow to reach regimes well inside the ordered 
phase, as shown, e.g., by the variation with F of \I/ m , the 
most probable value of the order parameter (Fig. [8]i). In 
both experiments Ea and Eb, \I/ m reaches values of the 
order of 0.5 for F = 2.8, but starts falling back for smaller 



T values. In the experiment on the large system Eb, the 
rise a \I/ m seems steeper than in experiment Ea, but the 
two curves are in fact difficult to compare given that 
the two experiments possess different averaged packing 
fractions in their ROI. 

Thus it is not clear from the variation of the order 
parameter alone whether our system, in its most ordered 
regime (F = 2.8), has reached the ordered phase or is still 
influenced by the transitional region. A direct measure- 
ment of the spatial and temporal correlation functions 
of the Eulerian velocity field suggests, though, that the 
ordered phase is indeed reached. The Eulerian veloc- 
ity field is here defined on a grid with a square mesh of 
size lxl. At each time step the local Eulerian velocity 
is the average of the instantaneous velocity of the par- 
ticles, the center of which is inside one element of this 
mesh. For small-enough F values, correlations in space 
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FIG. 9. (color online) Finite-size effects, (a) Spatial correlation C v (r) of the Eulerian velocity field v(f,t) for F = 2.8 for the 
small (black curves) and the large (red curves) systems at three different vibration amplitudes F as indicated in the legend, (b) 
Temporal correlation C v (t) of v(r,t) for F = 2.8 for the small (black curves) and large (red curves) systems at three different 
vibration amplitude F as indicated in the legend, (c) Fluctuations AiV 2 of the number of particles as a function of the average 
number of particles N in boxes of increasing size for the small (black o) and large (red +) systems. 



remain strong up to scales where the correlation function 
is cut-off by the system size, signaling that order is estab- 
lished on scales as large as possible (Fig. [9^l). Similarly, 
the temporal autocorrelation function decays very slowly 
(logarithmically?) for small F values and is only cut-off 
at very large timescales (Fig. (9b). 



D. Giant number fluctuations 

The orientationally-ordered phases of active matter 
systems are endowed with specific yet universal proper- 
ties such as long-range correlations, superdiffusion, and 
strong density fluctuations [TJ[T8]. The latter, which go 
under the vocable "giant number fluctuations", are the 
better known, probably because they are rather easy to 
measure: one simply records, in subsystems containing 
on average n particles, the variance An 2 of the fluc- 
tuations in time of this number. For normal, equilib- 
rium, systems, and/or for non-interacting particles An 2 
scales like n at large n (as long as, of course, n <C 
iV, the total number of particles in the system). But 
for orientationally-ordered active systems, one expects 
An 2 ~ n 1 with 7 > 1. 

A generic simple argument by Ramaswamy et al. [47] , 
links a q~ 2 divergence in the structure factor, calcu- 
lated in some linear, mean-field approximation, to 7 = 
1/2 + 1/d where d is the space dimension. This is ex- 
pected to hold for so-called "active nematics", and was 
indeed measured there [21] , but not for collections of po- 
lar particles as considered here. In this case, fluctuations 
must be taken into account, which has been done in the 
RG calculation by Toner and Tu [27 . This calculation, 
which relies on assumptions [28], predicts 7 = 8/5 for 
d = 2, a number consistent with numerical results ob- 
tained on the Vicsek model 22 . 



We measured number fluctuations in both experiments 
Ea and Eb for F = 2.8. Given the variations of local aver- 
age packing fraction described above, we used for n the 
mean number of particles actually recorded in each box 
in which particles were counted. This actually allowed 
us to use boxes slightly outside the ROI for better statis- 
tics without any detectable problem. As expected, An 2 
varies algebraically with n until it levels off at large n due 
to finite-size effects (Fig. [9^). The effective scaling expo- 
nent, measured over two decades, is 7 = 1.40 ± 0.01 for 
experiment Ea and 7 = 1.44 ± 0.01 for experiment Eb. 
This is close and below the expected theoretical value 
1.6, in agreement with the observation that the asymp- 
totic value is approached from below as the system size 
increases 21] . 



E. Intermittent confinement-induced milling 

Experiment Ec illustrates some of the consequences 
of increasing the density of particles. With a nominal 
packing fraction = 0.58, the ROI being fixed as previ- 
ously to a diameter of 20, we observe an effective pack- 
ing fraction inside the ROI </>roi = 0.19 suggesting a 
stronger concentration of the particles along the bound- 
aries, which cannot be "ignored" anymore: the time se- 
ries of \I/ in the ordered state observed at F = 3 show an 
intermittent behavior between values of the order of 0.6 
and values in the range [0,0.3] (Fig. [l0^-b). These low 
values do not indicate a disordered state. Rather, they 
are the trace of a milling configuration, i.e. a system- 
wide vortex obviously resulting from the domain shape 
(Fig. 10 >e and Supplementary movie [48 ). This can be 
quantified by calculating the macroscopic angular mo- 
mentum Fl(t) = (ni(tje0.(t))i, where e^(t) is the az- 
imuthal unit vector associated to 0$, the angular position 
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FIG. 10. (color online) Emergence of confinement- induced milling and its relation to polar order, (a) Temporal evolution of 
the ortho-radial ordering parameter IT(t) and the polar ordering parameter ^(t). (b) Zoom on a time window during which 
two inversion processes occur (vertical dashed lines) (c,d,e) Successive snapshots of the system with, from left to right, the 
inversion process from a clockwise ((II) < ) to an anti-clockwise (II) > large scale vortex. The packing fraction cj) = 0.58. 



of particle i in polar coordinates. Time series of \I/ and 
II reveal that the system switches intermittently between 
clockwise and counterclockwise milling: II changes sign 
from time to time, but never dwells long around zero. 
The usual order parameter in fact, takes large values 
mostly during these reversal events (Fig. flOk-b). 



The transition scenario we have observed in experi- 
ments Ea and Eb and described in the previous section 
is thus strongly modified. Here when T is decreased, 
one first observed the onset of collective motion, as in 
the previous case, but soon they organize in system size 
vortices and the transition towards polar order is inter- 
rupted. Further investigations are required to investigate 
the interplay between these two transitions. But one can 
already conclude that confined active flows are like any 
confined flows: boundaries cannot be ignored, and it is 
hard to disentangle the intrinsic "bulk" properties of ac- 
tive fluids from those resulting from the inevitable pres- 
ence of boundaries. 



VI. SUMMARY AND DISCUSSION 

Our results constitute, to our knowledge, the first 
study of a well-controlled experimental system in which 
"self-propelled" objects are able to move collectively over 
scales as large as the system size. The dynamics of 
our particles, although nominally three-dimensional and 
complicated, is well accounted for by a two-dimensional 
description in terms of persistent random walks. Their 
binary collisions are not simple inelastic aligning ones, 
but are spatiotemporally extended events during which 
multiple actual collisions happen, leading eventually to a 
weak but clear effective alignment. 

At the collective level, we showed how to avoid the 
accumulation of particles near the boundary walls by 
adopting flower-shaped arenas whose petals help reinject 
particles in the bulk. This trick is of course not perfect, 
and we do observe higher average densities near the walls, 
so that take must be taken in defining a region of interest 
in which the density is nearly constant although differ- 
ent from the nominal packing fraction. We presented 
the results of three sets of experiments conducted near 
"optimal" densities such that well-developed collective 
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motion is observed, in which orientational order is estab- 
lished on scales comparable to the system size. In these 
most ordered regimes, we recorded clear, unambiguous 
evidence of "giant number fluctuations" , a signature fea- 
ture of orient ationally-ordered phases in active matter. 

From the wealth of numerical studies of models of ac- 
tive matter, where periodic boundary conditions are of- 
ten adopted, it is easy to forget about the inevitable 
role played by walls when investigating collective mo- 
tion. Here, in addition to the subtle effects about the 
effective packing fraction recorded in the region of inter- 
est, we showed that at densities slightly larger than those 
used in the regimes where giant number fluctuations were 
recorded, the boundaries "come back" in the problem by 
driving macroscopic vortical flows. 

The experiments we conducted have also shown some 
clear limitations of our setup: Even though we used two 
domain sizes, finite-size effects remain strong and diffi- 
cult to estimate, which rules out any statement about 
the (asymptotic) nature of the transition to collective 
motion. Also, because our particles cease to move below 
some vibration amplitude value, we could not explore 
regimes "deep in the ordered phase" . Thus, to be true, it 
is hard to disentangle, in the strong density fluctuations 
measured, those intrinsic to the ordered phase from those 
linked to the proximity of the transition (not to men- 
tion again finite-size effects). Moreover, recent theoreti- 
cal progress has shown that the ordered phase, in a pa- 
rameter region bordered by the transition, is generically 
endowed with the emergence of high-density high-order 
bands [22j|24]. Here, we did not record any evidence of 



such bands, but one can argue that this is only because 
our system is too small. If these bands exist, then there 
is no theoretical reason to believe that the fluctuations 
exponent 7 = 8/5 predicted by Toner and Tu should 
be found, as their result applies to the (fluctuating) ho- 
mogeneous phase typically observed deep enough on the 
ordered size of the transition, but not when bands are 
present. 

The above remarks stress the difficulties inherent to 
systems of vibrated granular particles. In order to ap- 
proach asymptotic regimes closer, one should use larger 
domains and/or smaller particles, but both are techni- 
cally difficult. Reverting to the study of large swarms 
of small robots [49], or to biological motility assays [5], 
would suppress limitations in domain size, but at the 
price of some lesser control. Another approach, that we 
are currently pursuing [50], is to study in silico the sys- 
tem sizes and boundary conditions experimentally inac- 
cessible by constructing a numerical model faithful to 
every experimental finding. Obviously, such a model can- 
not be yet another Vicsek-style model. It has to incorpo- 
rate ingredients such as the treatment of the polarity vs. 
velocity dynamics, and to account for the complex series 
of collisions typically observed when two particles meet. 
It is our hope that such a model could give us access to 
the truly large-scale behavior of our system, and maybe, 
unveil new asymptotic properties. 

We thank M. van Hecke for advice in the design of our 
system, V. Padilla and C. Gasquet for technical assis- 
tance, and E. Bertin for enlightening discussions. This 
work was supported by the French ANR project Dy- 
CoAct. 
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